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Abstract —It is a great challenge to evaluate the network 
performance of cellular mobile communication systems. In this 
paper, we propose new spatial spectrum and energy efficiency 
models for Poisson-Voronoi tessellation (PVT) random cellular 
networks. To evaluate the user access the network, a Markov 
chain based wireless channel access model is first proposed 
for PVT random cellular networks. On that basis, the outage 
probability and blocking probability of PVT random cellular 
networks are derived, which can be computed numerically. 
Furthermore, taking into account the call arrival rate, the path 
loss exponent and the base station (BS) density in random 
cellular networks, spatial spectrum and energy efficiency 
models are proposed and analyzed for PVT random cellular 
networks. Numerical simulations are conducted to evaluate 
the network spectrum and energy efficiency in PVT random 
cellular networks. 

Index Terms —Random cellular networks, Poisson-Voronoi 
tessellation, Markov chain, spectrum efficiency, energy effi¬ 
ciency 


I. Introduction 

W ITH the development of the new generation cellular 
networks, more than 2,000,000 base stations (BSs) 
have been deployed until 2012. Report from the US shows 
that the energy consumption of an on-grid BS costs 3,000 
dollars per year on the average while the average energy bill 
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for running an off-grid BS reaches 30,000 dollars per year 
|H1 . Considering above energy cost and the carbon footprint 
in telecommunication systems, energy efficiency of cellular 
networks has received much attention from service providers 
G) . Therefore, it is important to investigate the spectrum and 
energy efficiency of cellular mobile communication systems 
0 - 0 . 

In cellular networks, a significant amount of power is 
consumed for transmitting signal over wireless channels. 
Therefore, the wireless channel model is a key metric for the 
performance evaluation of cellular networks. Gilbert-Elliott 
channel models were recommended as the most typical 
discrete channel model, where channel states are divided 
into the good or bad state 0. Based on the Gilbert-Elliott 
channel model, the network throughput and the bit error 
rate were investigated for wireless networks with bursty and 
noisy channels 0. Furthermore, the delay and throughput 
of cellular networks were analyzed based on Gilbert-Elliott 
channel models 0 ,0. Recently, finite state Markov mod¬ 
eling of fading channel was extended from wired network 
to wireless networks and used for wireless network per¬ 
formance analysis ED. Particularly, the received signal-to- 
noise ratio (SNR) was partitioned into a finite number of 
states and then a finite state Markov chain (FSMC) was 
constructed for modeling Rayleigh fading channels fill . 
Both numerical and simulation results have shown that the 
FSMC channel model provides an accurate model for the 
real channels. Based on a three-dimensional Markov chain 
model, a pico-cellular airport traffic model was proposed 
in oa which models Engset distributed fresh call arrival 
process and the general distributed handoff process with dy¬ 
namic channel allocation. The energy efficiency of fixed-rate 
transmissions from Markov sources over Rayleigh fading 
channels was investigated in na. On that basis, a closed- 
form expression of the minimum energy per bit was derived. 
Using the Markov chain to model the channel access in 
femtocell networks, energy and spectrum efficiency models 
were proposed in M for a two-tier femtocell network with 
partially open channels. 

Although one-dimensional linear model and two- 
dimensional lattices model including square lattices and tri¬ 
angular lattices were widely used in the modeling of cellular 
networks Ha, the assumption of regular deployment of 
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BSs in a plane is not consistent with real BSs deployment. 
Ignoring these structural fluctuations of the geometric ob¬ 
jects may cause significant bias on the evaluation of the 
cellular network performance lfl6l . To evaluate the geomet¬ 
ric structure of BSs coverage, the Poisson point process 
has been presented for modeling of BSs spatial structure 
03, Cl- Using realistic BSs locations collected from the 
Ofcom region UK, the Poisson point process has been 
justified for modeling of BS locations in cellular networks 
m. Moreover, in many cellular network scenarios only 
a statistical description of BSs locations is available. For 
these scenarios, the stochastic geometric model provides a 
promising solution for modeling and performance analysis 
of cellular networks. Particularly, in ED, Win et al. intro¬ 
duced a mathematical framework for the characterization of 
network interference in which the interferers are scattered 
according to a spatial Poisson point process and are oper¬ 
ating asynchronously in a wireless environment subject to 
path loss, shadowing and multipath fading. Using physically 
realistic stochastic models for BS’s locations in multi-tier 
heterogeneous networks, a general signal-to-interference- 
and-noise ratio (SINR) distribution was derived in [20] for 
calculating the probability of the user being able to camp on 
a macrocell or an open-access femtocell in the presence of 
closed subscribe group femtocells. Assuming mobile users 
and BSs in different tiers of heterogeneous cellular networks 
are independent Poisson point processes, a tractable model 
of downlink heterogeneous cellular network consisting of 
K tiers of randomly located BSs was developed in CTR . 
By assuming spatial Poisson distribution of transmitters in 
the primary network, the achievable transmission capacity 
of a secondary cognitive mesh network sharing the uplink 
spectrums with cellular users under the outage probability 
constraints of both the primary and secondary systems was 
investigated in l22l . Using a stochastic geometry based 
model and different sleeping policies, the success proba¬ 
bility and energy efficiency in homogeneous macrocell and 
heterogeneous K-tier wireless network were derived in [23]. 
In l24l . based on the spatial Poisson point process, the 
distribution of SINRs and mean achievable rates of both 
non-cooperative users and cooperative users were derived 
to evaluate the energy efficiency of relay-assisted cellular 
networks. 

In all aforementioned studies, Markov chain models and 
stochastic geometry tools are successfully used for inves¬ 
tigating cellular networks, respectively. However, the study 
combining the Markov chain model with random cellular 
networks is surprisingly lacking in the open literatures. We 
started an initial study and published the initial result in 
(25). One of the major obstacles in combining the Markov 
chain model with random cellular networks is how to model 
the user access in stochastic geometry cellular scenarios by 
Markov chains. To overcome this obstacle, we propose a 


Markov chain based channel access model for a Poisson- 
Voronoi tessellation (PVT) random cellular scenario. Fur¬ 
thermore, the Markov chain based channel access model for 
one PVT cell has been extended to the whole PVT random 
cellular network using the Palm theory. As a consequence, 
the spatial spectrum and energy efficiencies of PVT random 
cellular networks has been analyzed in this paper. The 
contributions and novelties of this paper are summarized 
as follows. 

1) Although the Markov chain model and the stochastic 
geometry theory have been widely used in cellular 
networks, it is still a great challenge to integrate the 
Markov chain model into the analysis of stochastic ge¬ 
ometry cellular networks. Based on the PVT cellular 
scenario, a Markov chain based channel access model 
is first proposed and integrated into random cellular 
networks. 

2) From the proposed wireless channel access model, 
we derive the exact outage probability and blocking 
probability for PVT random cellular networks with a 
spatial Poisson distribution of interfering transmitters, 
taking into account fading and shadowing effects over 
wireless channels. 

3) Using the outage probability and blocking probability, 
a new spatial spectrum and energy efficiency model 
of PVT random cellular networks is developed for 
performance evaluation. 

4) We study the spatial spectrum and energy efficiency of 
PVT random cellular networks in details and present 
some interesting observations. 

The remainder of this paper is outlined as follows. Sec¬ 
tion [III describes the system model. In Section HIH a wireless 
channel access model based on Markov chain is proposed 
and on that basis the outage probability and blocking 
probability of PVT random cellular networks are derived. 
Furthermore, the spatial spectrum and energy efficiency 
models of PVT random cellular networks are presented and 
evaluated in Section UV] Finally, Section [V| concludes this 
paper. 

II. System Model 

Assume that both mobile users (MUs) and BSs are 
located randomly in the infinite plane M 2 . Moreover, 
locations of MUs and BSs are modeled as two inde¬ 
pendent Poisson point processes, which are denoted as 
Qu={xi :i = 0,1,2,...} and 0 B = {yj,j = 0,1,2,. 
where Xi and yj are two-dimensional Cartesian coordinates, 
denoting locations of the ith MS MUi and the jth BS 
BSj, respectively. The corresponding intensities of the two 
Poisson point processes are A u and A#, respectively. For 
a cell, the interference among MUs is ignored and the co¬ 
channel interference is assumed to be mainly contributed by 
adjacent BSs. 
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A. Wireless Propagation Environments 

For an MU xi G ©u and the associated BS yj G 0#, 
the channel gain L y .(xi), defined as the ratio between the 
received power P r _ Xi at an MU and the corresponding 
transmission power P y . from its associated BS, is 


L Vl {xi) 



K-S yj (xj) 

L(\\Vj—Xi{\) \\yj - Xi\\ b 


where K is a constant depending on antenna gains, 
SyAxi)=Uk Z t( x i) is the impact of fading and shadow¬ 


ing effects on the received signal power in wireless prop¬ 


agation environments; 


s; | Zy. (pci) k G Z + | are independent 
random variables (RVs), which account for propagation 
effects such as fading and shadowing, where z+ is the 
positive integer set, the values of k and the distribution 
of ^Zy J . (^i)j are general enough to account for various 
propagation scenarios, including the following: 


1) Sy^X^ = Zy.(Xi) = 1: without fading and shadow¬ 
ing effects; 

2) Zy (x^ follows an exponential distribution: Rayleigh 
fading propagation scenario; 

3) S Vj (xi) = Zy.(xi) = ? 2 with ? 2 - r(m, ^), 
where V (m, —) is a Gamma distribution with 1 mean 
and variance of Nakagami-m fading propagation 
scenario; 

4) Sy.[xi) = Zy. (xi) = e 2aG with G ~ 
Gaussian (0,1), where a is the shadowing coeffi¬ 
cient: log-normal shadowing propagation scenario; 

5 ) S yj (xi) = n k Zy. (xi) = Z^(xi)Z^(xi) = c 2 -e 2aG : 
Nakagami-m fading and log-normal shadowing prop¬ 
agation scenario. 


L(\\yj — Xi ||) accounts for the far-field path loss, which is 
a function of the distance between the MU located at Xi 
and the BS located at yi , denoted as \\yj — Xi\\ b , where b is 
the path loss exponent. In general, the path loss exponent is 
environment-dependent and can approximately range from 
0.8 (e.g., hallways inside buildings) to 6 (e.g., buildings 
without line of sight) (T9), ll26l . 


B. User Association Scheme 

In realistic cellular networks, the MU receives different 
pilot signals from adjacent BSs and then associates with a 
BS based on a specified user association scheme. In recent 
studies, different user association schemes were adopted 
for evaluating the performance of cellular networks, for 
example, the nearest BS association scheme was adopted 
in GQ), G 3 , Gz), the highest SINR association scheme 
was adopted in 1211 . |[28lL f29lL and the maximum received 
signal power association scheme was adopted in lf30l — lf32l . 
When an MU located at Xi is associated with a BS located 
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Fig. 1. An illustration of the PVT random cellular network: blue 
points are BSs and red points are MUs. 


at yj, a general association scheme is expressed as [33, 

GD- 

y* = arg max T Vj \\yj - x t || _b , (2) 

VjE'&B 

where T yj is the weight association metric, which is de¬ 
pendent on the particular user association scheme. In this 
paper, the MU is assumed to be associated with the nearest 
BS in the plane. In this case, the weight association metric 
is configured as T yj = 1. Accordingly, (2) can be rewritten 
as 

y* = arg max \\yj - Xi\\~ b . (3) 

yjee B 


Based on the PVT method, the coverage of BS located at 
yj is defined as [34 ] 


C Vj = (2/ e K 2 : II y - yj\\ < ||y - y;|| ,Vy; e 0s and 


Vi 7 L Vj\- 


(4) 


Based on ©, an illustration of the PVT random cellular 
network is depicted in Fig. Q] where each cell is denoted 
as C y . . In Fig. [I] the cells are rounded by blue lines, the 
BSs are denoted by blue points and the MUs are denoted by 
red points. Such stochastic and irregular-topology network 
forms the so-called PVT random cellular network. 

Despite of its apparent complexity, an outstanding prop¬ 
erty of PVT random cellular networks is that the geometric 
characteristics of any cell C Vj coincide with that of a typical 
PVT cell C or i where the corresponding BS is located at a 
fixed position, e.g., the origin, according to the Palm theory 
EH, ( 22 . This feature implies that the analytical results 
obtained for a typical PVT cell C or i can be extended to 
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the whole PVT random cellular network. Without loss of 
generality, the coverage of the typical cell C or i is expressed 
as 

Cori = {y e R 2 : || V - ori\\ < ||y - yi\l ,\/yi G © B and 
Vi ^ ori}, 

( 5 ) 

where ori = (0,0) is the origin point in the plane R 2 . 

C. Wireless Channel Allocation Strategy 

In this paper, a centralized channel allocation strategy is 
used for PVT random cellular networks. The traffic channels 
are allocated for a MU after this MU is associated with a 
specific BS. Furthermore, every MU is assumed to have 
the perfect channel state information (CSI) of all channels 
available to its BS. This means the associated MU can 
get the SINR over all channels. When the SINR value 
of a channel is larger than or equal to the threshold 70, 
this channel is available for the associated MU and will 
be marked as 1. When the SINR value of a channel is 
less than the threshold 70, this channel is unavailable for 
the associated MU and will be marked as 0. All available 
channel states will be feedbacked to the control centres of 
PVT random cellular networks, and then the control centres 
will allocate BS channels to all associated MUs based on 
the set of available channels. 

III. Models of PVT Random Cellular Networks 
A. Markov Chain Model of PVT Random Cellular Networks 

Without loss of generality, a typical cell C or i in the 
PVT cellular network is selected for Markov chain based 
modeling and performance analysis in this paper. The arrival 
rate of calls to the system is assumed to follow a Poisson 
distribution with mean A. The user session duration Ts and 
the cell dwelling time Td are assumed to be governed by ex¬ 
ponential distributions with mean /i and l/Tp respectively, 
where Td is the mean cell dwelling time. The channel 
holding time Th is the minimum of the user session duration 
and the cell dwelling time, i.e., Th = min (T^, Td), with 
mean 7 = fi + 1/Td . Unlike the discrete Gilbert-Elliott 
channel model used in many articles 0 , 0 , in this paper 
the continuous Gilbert-Elliott channel model is used to 
model the transitions between the available channels and the 
unavailable channels. For traditional Markov chain models 
of wireless channel access in cellular networks, the state 
change of unavailable BS channels has been modeled by a 
Markov chain and the state dwelling period of unavailable 
BS channels is governed by an exponent distribution 0 , 
0 . Considering that the number of available BS channel 
is determined by the total BS channel number minus the 
unavailable BS channel number in a cell, the state dwelling 
period of available BS channels is also assumed to be 
governed by an exponent distribution. The transition rate 
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Fig. 2. Markov chain transition diagram. 

from the unavailable channel to the available channel is 
denoted as a and the transition rate from the available 
channel to the unavailable channel is denoted as /3. The 
two parameters will be determined later. 

A two-tuple of nonnegative integers (m,n) is used to 
denote the network states, where m is the number of 
occupied channels and n is the total number of channels 
available to be allocated for MUs in the typical cell C or i , 
including the channels that have already been allocated (i.e. 
occupied). Let C be the maximum available channels in the 
typical cell C or i. Obviously, the constraint m < n < C 
has to be observed. The Markov Chain transition diagram 
is illustrated in Fig. O 

The Markov Chain state transitions in Fig. [2| are described 
as follows: 

1) (m,n) -» (m + l,n): When a new call arrives, the 
number of occupied channels is increased by 1 if the 
number of occupied channels is less than the number 
of available channels, i.e., m < n. 

2) (m, n) -A (m — 1, n): When a call has been success¬ 
fully serviced, the number of occupied channels is 
reduced by 1 if this occupied channel is released. 

3) (m,n) -A (m, n + 1): An unavailable channel be¬ 
comes available to be allocated for mobile users due 
to the time-varying interference, thus the number of 
available channels is increased by 1 . 

4) (ra, n) -A (m,n — 1): An available channel becomes 
unavailable to be allocated for mobile users due to 
the time-varying interference, thus the number of 
available channels is reduced by 1 . 

Based on the Kolmogorovs criteria 03, a Markov pro¬ 
cess is reversible and its stationary state distribution exists 
if its transition rates satisfy 

A (Si, S 2 ) • A (S 2 , S 3 ) •... • A (S u - lt S u ) • A (S u , S x ) 

= A (Si, S u ) • A (S u , S u -i) •... • A (S 3 , S 2 ) • A (S 2 , Si) 

( 6 ) 

for any finite sequence of states Si, S 2 , • • • , S u GS, where 
E is the set of all possible states in this Markov process 
and A (S^, S&) denotes the transition rate from state Sh to 
state Sfc. It is observed that states in Fig. [2| meet © and the 
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" Col • 7 r( 0 , 0 ) = /3 • 7 r( 0 , 1 ) 

[/3 + A + (C — n)\ • 7r(0, n) = (C — n + l)a • 7r(0, n — 1) + 77 • 7r(l, 72 ) + (72 + l)/?7r(0, n + 1) 

/or 0 < n < C 

[C/3 + A] • 7 r( 0 , C) = a • 7 r( 0 , C — 1 ) + 7 • 7 r( 0 , C) 

[C/3 + A + 77277 ] • 7r(l, C) = a • 7r(m, C — 1) + (m + 1)77 • 7r (772 + 1, C) + A • 7r(m — 1, C) 

< for 0 < m < C ( 7 ) 

Cr\ • 7r(C,C) = A-tt(C-1,C) 

[(C — rri)a + 77277 ] ' 2 / 2 ) = A • it{m — 1, m) + (rri + l)/3 • 7 t(t72, m + 1) 

/or 0 < 772 = n < C 

[n /3 + (C — 72)<Y + A + 77277 ] ■ 7r(^j 72 ) = (C + 1 — 72 ) O' • 7r(772, 72 — 1) + (72 + l)/3 • 7r(772, 72 + 1) + 

(t 72 + 1)77 • 7r(772 + 1, 72) + A • 7r(772 — 1, 72) /or 0 < 772 < 72 < C 


stationary state distribution exists as a result. For an easy 
illustration, four states illustrated with purple dash line are 
taken as an example. State transition rates on the left and 
right side of © are placed beside the red arrows and blue 
arrows, respectively. In this way, the product term on the 
left side of © is a • 77 • C/3 • A and 77 • a • A • C/3 on the right 
side. Obviously, <Y-? 7 -C/ 3 -A = ? 7 -auA-C/ 3 . Similarly, we 
can verify other finite sequences of states and find that © 
holds. According to the Markov chain state transitions in 
Fig. |2l global equilibrium equations are listed as ©. Based 
on the queuing theory, the stationary state probabilities are 
derived as follows 


7r (772,72) 




( 8 ) 


where (^) is the binomial coefficients meaning the number 
of ways of picking n unordered outcomes from C possibil¬ 
ities. 

Based on the Gilbert-Elliott channel model ( 8 ), (9J, the 
probability of a channel being unavailable e is defined as 


£ = 


P 

+ /3’ 


(9) 


The ratio between a and /3 needs to be calculated in ®. 
Based on a//3 equals to (1 — e) /e. In this paper, the 
channel is unavailable when the SINR value of channel is 
less than the given threshold. In this case, the probability 
of unavailable channel is equal to the expected value of the 
outage probability with different number of interferers. e 
then is expressed as follows 


£= J im 'YhPont (7°, A) -p) Nl A , (10) 

AT/—>-00 z ' \ A / 

A = 1 x 7 

where 3Vj is the maximum number of interferers in 0 ^; 
A(A < Ni) is the possible number of interferers aggregated 


at a wireless channel; p out (70, A) is the outage probability 
conditioned on that there are A interferers; p is the proba¬ 
bility that a channel is being used at present, i.e., 

n C 

(11) 

m =0 n =0 

It is observed that the stationary state probability a//3, the 
probability of a channel being unavailable 5 and the prob¬ 
ability of a channel being used p are mutually dependent 
and correlated with each other in ®, ©, TO and CD. 
Fortunately, above probability values can be obtained by 
solving linear equations if we know the outage probability 
Pout ( 70 , A). In the next paragraph, we compute the value 
of the outage probability. 


B. Outage Probability Model 

Considering that in the typical cell C or i , the MU located 
at Xi is associated with the BS located at 77 and co-channel 
interference from adjacent BSs exists. Based on the general 
wireless propagation environments described in Section [III 
the SINR at the MU located at Xi is expressed as follows 


SINR Vj (x^ 


Kys^HL ii )}- 1 

cr2 + I £ 


( 12 a) 


with 


7 = E K' ■ S m (xi){L (\\ yi ~ XiW)}- 1 , 
yiee B \{yj}, 

( 12 b) 

where is the aggregate interference received at the MU 
conditioned on that there are A interferers in total; #{•} 
is an operation which counts the number of elements in a 
set; cr 2 is the Gaussian noise power. Assume that all BSs 
transmit with the same power and as a result, K' = KP y ., 
Vj € $£• 
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Because the sum of the outage probability p out (70, A) 
and the success probability p suc ( 70, A) equals to 1, 
Pout (70, A) is expressed by 

Pout ( 7 O 5 A) = 1 Psuc ( 70 ? A) . (13) 


When the MU is assumed to be associated with the 
nearest BS, the success probability of a MU located at 
E C or i is derived as (fl4l ) . where E{*} is an expectation 
operation, g (t) is the probability distribution function (PDF) 
of (cr 2 +/^) and 1 {•} is an indicator function, which 
equals to 1 when the condition inside the bracket is satisfied 
and 0 otherwise. Moreover, the success probability can be 
further derived by using Parseval theorem |[T6l given below. 

Parseval theorem : Denote the Laplace transforms 
of two square-integrable functions / (t) and g (t) as F (lj) 
and G (lj), respectively, i.e., f (t) = A_ p (cc) e luJt dw 
and g (t) = A_ G (u) e luJt duj, then they have the 
relationship as follows 


+00 


1 n°° 

f(t)g(t)-dt=—J F (cu) G (cu) ■ doj. (15) 

Let the PDF of S y .(xi) be f (y) and assume that 
f (y) is independent of g(y). Based on (IT5l) and the 
assumption that the noise and interference are indepen¬ 
dent of each other, (m can be further derived as m 
where C a 2 (•) and Cj a (•) are the Laplace transforms 
of noise and interference /^, respectively. Assume that 

dwdy is 
tainea as 

Further with the assumption of dTSb . (fTTb can be further 
derived as 01, where Cs ( Xi ) (') I s the Laplace transform 
-u)K' • {70L(r)} -1 and (fl 9 l) is 


fo°° f-oo { ex P [ 2 ^wj/^{i(r)} - 1} 

finite, then the success probability can be obtainec 


of Su 


j yj (xi). Let s = 

expressed by (l20l) . 

The Laplace transform of the noise is expressed as m 


C a 2 (s) = e~ 


( 21 ) 


According to the definition of /pj in (fT2bl we follow 
a basic technique in [[38) to obtain the outage probability 
condition on having A interferers, which consists of two 
steps: 

1) Consider a finite network, say on disk of radius 
a centered at the origin, and condition on having 
a constant number of nodes in this finite area, for 
example A nodes. Assume that the nodes locations 
are independent identical distribution. 

2) Let the disk radius go to infinity, while keeping the 
node density, i.e. the ratio of the number of nodes to 
the network area, constant. 

Stepl: Conditioning on having A nodes (i.e., interferers) 
in the disk of radius a, the Laplace transform of aggregate 
interference is denoted by 


,a (s) = E { e~ sI ^ | #£? ( ori, a) = A} , (22) 


where B (ori, a) is the set of nodes located on a disk with 
radius a centered at the origin ori and j^B (ori, a) = A 
means that the number of nodes in B (ori, a) is A. These 
nodes (i.e., interferers) are uniformly distributed on the disk 
with radial density 


■MO 


2 r 
a 2 

0 


if 0 ^ r ^ a 
otherwise 


(23) 


Therefore, the Laplace transform jCja a (s ) is the product 
of the A individual Laplace transform, which is given by 
©. 

Assume that the path loss law is L (r) = r b , then the 
individual Laplace transform is derived as follows 

J Q ex P (—sK 1 ■ S yi {xi){L (r)} -1 ) dr 


a 

r a p r 

= — exp (- sK' ■ S Vl ( Xi)r~ b ) dr 

J 0 a 
2 f°° 

= - 2 / y ~ 3 ex P (~ sK> • Syi ( x i)v b ) dy 

a J x 

a 

= \[sK' • s yi { Xi )Y ■ r (-IsK'^SyAXi) 

where F (•) denotes Gamma function. 


(25) 


Step2: When a goes to infinity, (l25t is further derived 
by 

n CL c\ 

lim / — exp (-siF' -S Vl (xi){L dr 

a—>-oo Jq \ J 

= l -[sK’-Sy^Y ■ r (-|,o) 

= l -[sK’ ■ Sy^Xi)] 1 -T (--\ 

(26) 

Then Ci x . A (s) is is transformed into (l27l) . By sub¬ 
stituting df and d27l> into ©, the outage probability 
conditioned on having A interferers is derived by (l28l) . 


C. Blocking Probability with Rayleigh Fading 

When wireless signals are assumed to suffer Rayleigh 
fadings, the PDF of signal is expressed as 

f(x) = x/a 2 • exp (—x 2 /2ct 2 ) (0 ^ x ^ oc). (29) 

And the PDF of the signal power is expressed as 

f(y) = l/2cr 2 • exp (-y/2cr 2 ) (0 < y < oo). (30) 

It is obvious that the PDF of the signal power follows 
an exponential distribution. Without loss of generality, let 
Syi(xi) ~ exp(l). Then the Laplace transform of S yi (xi) 
is derived by 

^s Vl (xi) ( s ) = s + i’ 








IEEE TRANS. ON COMMUNICATIONS, VOL. XX, NO. Y, MONTH 2015 


7 


= P( 

SINR Vj (x^ > 7 o| Xi 

— p 1 

'( K'■ S Vj {xi){L (\\ yj - 

\ 

{ <T 2 +I Xt ,A 

=p( 

[s Vj (x i )>^ J .L(\\y : 

=H 

!'(-«<?• 

H 

[ 2 ttA B re _7rABr2 d 
l *7 r >0 


l -1 


> 7o 


Xi G C 0 


7 *r?" 


G 2 + 4l) 

r / \— 1 


37 G 




/*+00 


(14) 


p suc ( 70 , A) 

E <( [ 2ttA Bre-^^ 2 dr < 


/ r>0 


j_ C 00 

2tt J_ oq 


C a 2 (—2njuj) ■ £/a (— 2t: ju>) 


j ■■ 

J r> 0 


27rA^re 


— 7tA r7^ 


dr 


/ +00 

Hv ) 

-00 


r +00 




£ ct 2 (-27rjw) • £/a (-2njw) 


^2irjuj^--Syj (xi){ L (r)} 1 _ ^ 

27 xjuo 

o 2nju -y-{L(r)}~ 


duj 


- 1 


47T 2 JCcJ 


duj 


(16) 


( 70 , A) 

/ 27tA #re 

J r 


-TrAer rfr 


/r>0 


/ +00 

/(y) 

-00 




_ 7*+oo 

-- cfi/ / £ a 2 (— 2itjuS) • £/a (—27t joo) duj 

J — c 


/ : 

<7 r>0 


27rA^re 


r+oo 


— 7tAb£^ 


dr ■ 


47T 2 jCo’ 

A°° /_ + ~ / (y) • * d y _ / (y) ^ 


' —OO 

1-1 f + oo 

-00 


47r 2 jci; 


dec 


/ f 


£+00 
I r>0 J —00 


C a 2 ( — 27Tj(jj) • £/A ( — ‘lltjuS) duj 

x i 

50 

/_ + 7 / (y) • - /_ + ~ / (y) dy 


47r 2 jcj 


I At 2 (—27Tjr'u;)| • £/a (~2njw) 


(17) 


dccdr < 00 . (18) 


and then 


E 


{(Sv,te)) ¥ } 


= r 1 


2 A 

IT 


(32) 


Moreover, the Laplace transform of the aggregate interfer¬ 
ence is derived by 


Ci t , a (s)= [- h (sK')- b -V 


into 


■r 1 


will be blocked due to a lack of sufficient channel resource. 
Therefore, the blocking probability is derived by 

( \ ST ^ f ^ ^ ^ ^ ^ a 

p b = ^ w(m,n) = ^ -(-) ^ 

m=n^C 


with 


2 A\ 
(33) 

, the outage probability is 


i=j^C 

a 1 — e 


X\dJ m - \ n J \P 


P 


Nj 


By substituting 
derived as (HQ). 

Assume that a call will be dropped if the call can not be 
served immediately in PVT random cellular networks. As a 
consequence, if the total number of active MUs exceeds the 
number of available channels in the typical cell C or i , a call 


= 77oo E p ° ut Go, a) (7 b A (i - p) 

A = 1 


Nl\ A/-. \N T -A 


(35a) 

(35b) 

(35c) 


Furthermore, the mean sojourn time of a MU can be 
derived by Little’s theorem m 
„ iV 
D= A’ 


(36) 
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c 

where TV = ^ m • 7r (ra, n) is the mean number of served 

m=l 

MUs in the typical cell C 0 h. 

D. Performance Analysis 

Assuming Rayleigh fading channels, the blocking prob¬ 
ability and the mean sojourn time of PVT random cellular 
networks can be numerically computed. Following the simu¬ 
lation configuration in 13, ED and m, default parameters 
for a PVT random cellular network are configured as 


follows: BS intensity is X B = 0.2 per square kilometers; 
the maximum number of available channels in a typical cell 
C or i is C — 20; the arrival rate of calls is A = \minute~ x \ 
the BS transmitting power is P yi = 30dBm; the mean noise 
power is cr 2 = OdBm; the path loss exponent is b = 4; 
the antenna gain is K = 31.54dB for an urban microcell 
environment. 

Fig. a shows the call blocking probability with respect to 
the SINR threshold considering different maximum number 
of channel numbers C in a typical cell C or i. When the 
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(35d) 



Fig. 3. Blocking probability versus SINR threshold with different 
maximum channel numbers C. 



Fig. 5. Mean sojourn time versus SINR threshold with different 
A. 


0.9 



1 5 10 15 20 25 30 35 40 45 50 


SINR threshold y 0 

Fig. 4. Blocking probability versus SINR threshold with different 
b. 


maximum number of available channels is fixed, it is 
observed that the blocking probability increases with the 
increase of SINR threshold. The reason is that the successful 
decoded signal is reduced for MUs when the SINR threshold 
is increased. When the SINR threshold is fixed, the blocking 
probability increases with the decrease of the maximum 
available channels number in the cell C or i • That is because 
a call is more likely to be dropped when channel resource 
is insufficient. 

Fig. [4] illustrates the blocking probability versus the SINR 


threshold for different path loss exponents b in the typical 
cell C or i • When the SINR threshold is fixed, the blocking 
probability decreases with the increase of the path loss 
exponent. It is well known that the path loss exponent affects 
both desired signals and interference signals. However, 
these curves imply that the path loss exponent has a more 
significant attenuation on the aggregated interference in 
PVT random cellular networks. 

Fig. [5] shows the mean sojourn time versus the SINR 
threshold for different arrival rates of calls. When the SINR 
threshold is fixed, the mean sojourn time decreases with the 
increase of arrival rate of calls. When the arrival rate of calls 
are fixed, the mean sojourn time decreases with the increase 
of SINR threshold. That is because the number of available 
channels decreases with the increase of SINR threshold. As 
a result, the service ability of cellular network is decreased 
and then the mean sojourn time of a MU staying in the PVT 
random cellular network is decreased. 

IV. Spatial Spectrum and Energy Efficiency of 
PVT Random Cellular Networks 

In this section, we further evaluate the spatial spectrum 
and energy efficiency of PVT random cellular networks. 

A. Spatial Spectrum and Energy Efficiency 

Assume that the bandwidth of a typical cell C or i is B , 
the throughput of the typical cell C or i is then given by 
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T- throughput — (1 Pb) R ' E |log 2 (l SINRy. (3?i))} * E m • 7r(m, n). (37) 

O^m^n^C 


Without loss of generality, the number of interferers on a 
wireless link of PVT random cellular networks is assumed 
as A. Therefore, the link capacity between a MU and the 
associated BS is defined as 

E {log 2 (1 + SINR Vj (x i ))} 

P (log 2 (1 + SINR yj (x,)) > t )dt 

f381 

P (SINR.y. (xij > 2 l - l)dt 

+oo 

[1 -Pout (2‘-l,A )]dt 

A simple proof of (l38l ) is given as follows. 

Proof: Consider a random value x which has continues 
PDF, then its expectation is given by 




E 


/ too 

xf{x)dx 

-oo 

/ 0 p-\~oo 

xf(x)dx + / xf(x)dx 
-oo J 0 


(39) 


Furthermore, the two terms of the right side of (l39l) are 
written as 


[ xf(x)dx = - f { f dy \ f{pc)dx 
7—oo 7—oo V 7 x J 

/ O py nO 

/ f(x)dxdy = - F(y)dy 
OO 7 — 00 7 — 00 

r+oo r+oo / rx \ 

J xf(x)dx = J (J dy) f(x)dx 

/> + 00 /* + 00 /* + 0O 5 

/ / f(x)dxdy = [1 - F(y)} dy 

Jo Jy Jo 


(40a) 


(40b) 


where F (•) denotes the cumulative probability function 
(CDF) operation. As we know SINR Vj (xi) ^ 0. Substitut¬ 
ing SINR yj (x^ into (140bk the result of (l38t is obtained. 
Based on (l37k the spatial spectrum efficiency of the PVT 
random cellular network is derived as <ED according to 
ll23l . where A b is the BS density of PVT random cellular 
networks. 

The energy efficiency of the typical cell C or i during the 
whole life time is derived as follows 


P = 


Rfotal 

Dtotal 


(42) 


where E tota i denotes the BS energy consumption in the 
life time and D to tai denotes the BS throughput in the life 
time. In the entire life time, the BS energy consumption 
includes the operation energy and the embodied energy ED. 


Moreover, the embodied energy includes the initial em¬ 
bodied energy consumed in factories and the maintenance 
embodied energy in the life time. Therefore, the BS energy 
consumption in the life time is expressed as 

Etotai = EfiMinit H“ E^Mma int H - EfiMopen (43) 

where EEMinit denotes the initial embodied energy, 
EEMma int denotes the maintenance embodied energy, and 
Eem oper is the operation energy, which is a linear function 
of the total transmission power over all occupied channels 
and is given as follows lt42l 


EsMoper = h ■ P c hl E m ' W ( m > n ) + k 
O^m^n^C 

X tlifetime •> 


(44) 


where P c hi is the transmission power over a wireless 
channel, tuf e time is the life time of a BS, h and k are 
linear coefficients of the total transmission power. The total 
throughput of BS in the life time is expressed as 


Etotai — tlif e time ’ 2 "throughput • (45) 

Therefore, the energy efficiency of the typical cell C or i is 
derived as (l46l) . Based on the Palm theory, the result of (f46b 
can be extended to the whole PVT random cellular network. 

B. Numerical Results and Discussion 

Based on the spatial spectrum and energy efficiency 
analysis, numerical results are illustrated in this subsection. 
Based on default parameters used in Section HIIlD, some 
parameters used for comparing grid cellular networks (i.e., 
regular hexagonal cellular networks) and PVT random cel¬ 
lular networks are configured as follows J35|, ATI . [42): 
the channel bandwidth is B = 0.1MHz, the embodied 
energy is configured as E E Minit + E E Mma int = 85GJ, the 
transmission power over a wireless channel is P c hi = lWatt, 
h = 7.84, k = 71.5, without loss of generality, the BS life 
time tufetime is configured as 1 year, e.g., 365 days. 

Fig. [6] illustrates the spatial spectrum efficiency of PVT 
and grid cellular networks with respect to the path loss 
exponent and the BS density in cellular networks, in which 
“PVT model” labels the PVT random cellular network 
and “Grid model” represents the regular hexagonal cellular 
network. When the BS density A# is fixed, numerical results 
of PVT and grid cellular networks consistently show that 
the spatial spectrum efficiency increases with the increase 
of the path loss exponent. When the path loss exponent b is 
fixed, numerical results of PVT and grid cellular networks 
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Fig. 6 . Spatial spectrum efficiency with respect to the path loss 
exponent b considering different BS densities As in PVT and grid 
cellular networks. 


consistently illustrate that the spatial spectrum efficiency 
increases with the increase of the BS density. Compared 
with PVT model results and grid model results in Fig. [6j it 
is shown that values corresponding to PVT random cellular 
networks are obviously less than values corresponding to 
grid cellular networks. This result is confirmed in m which 
showed that the average transmission rate of PPP random 
cellular networks is less than the average transmission rate 
of grid cellular networks. Considering that the PVT random 
cellular network forms a special case of the PPP random 
cellular networks, the spatial spectrum efficiency of PVT 
random cellular network is less than the spatial spectrum 
efficiency of grid cellular network when the transmission 
bandwidth is fixed in PVT and grid cellular networks. 

In Fig. [71 the effect of the call arrival rate A on the 
spatial spectrum efficiency of PVT and grid cellular net¬ 
works is investigated. When the path loss exponent is 
fixed, numerical results of PVT and grid cellular networks 
consistently demonstrate that the spatial spectrum efficiency 
decreases with the increase of the call arrival rate. When the 
call arrival rate increases, the blocking probability in the 
cell is correspondingly increased. The increase of blocking 
probability in the cell in turn reduces the spatial spectrum 



Fig. 7. Spatial spectrum efficiency with respect to the path loss 
exponent b considering different call arrival rates A in PVT and 
grid cellular networks. 



Fig. 8 . Spatial spectrum efficiency with respect to call arrival 
rate A considering different SINR thresholds 70 in PVT and grid 
cellular networks. 


efficiency of PVT and grid cellular networks. 

Fig. [8] shows the spatial spectrum efficiency with respect 
to the call arrival rate considering different SINR thresholds 
in PVT and grid cellular networks. When the call arrival rate 
in cellular scenarios is fixed, the spatial spectrum efficiency 
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Fig. 9. Energy efficiency with respect to call arrival rate A 
considering different SINR thresholds 70 in PVT and grid cellular 
networks. 


Fig. 10. Energy efficiency with respect to path loss exponent b 
considering different BS densities As in PVT and grid cellular 
networks. 


decreases with the increase of the SINR threshold. When the 
SINR threshold is fixed, there exist thresholds for different 
call arrival rates in cellular scenarios. Below the threshold, 
the spatial spectrum efficiency increases with the increase 
in the call arrival rate and above the threshold the spatial 
spectrum efficiency decreases with the increase in the call 
arrival rate. Numerical results of PVT and grid cellular 
networks consistently validate that there exist maximum 
spatial spectrum efficiency values considering different call 
arrival rates in cellular scenarios. 

The impact of the call arrival rate on the energy efficiency 
of PVT and grid random cellular networks is evaluated 
in Fig. 0 When the SINR threshold is fixed, there exist 
thresholds for different call arrive rates in cellular scenarios. 
Below the threshold, the energy efficiency increases and 
above the threshold the energy efficiency decreases with 
the increase in the call arrival rate. Numerical results of 
PVT and grid cellular networks consistently validate that 
there exist maximum energy efficiency values considering 
different call arrive rates in cellular scenarios. Therefore, to 
achieve an optimal spectrum and energy efficiency of cellu¬ 
lar networks, the call arrive rate and the SINR threshold in a 
cell should be considered carefully by telecommunications 
operators. 

Fig. [T0| depicts the energy efficiency with respect to the 
path loss exponent considering different BS densities in 
PVT and grid cellular networks. When the BS density of 
cellular networks is fixed, numerical results of PVT and grid 
cellular networks consistently confirm that the energy effi¬ 
ciency increases with the increase of the path loss exponent. 
With the increase of the path loss exponent, both the desired 
signal and the interference are exponentially attenuated over 
wireless channels. Since the distance between the interfering 
transmitters and the receiver is longer than the distance 
between the desired BS and the receiver, the interference 


will experience larger attenuation than the desired signal 
when the path loss exponent increases. Therefore, the outage 
probability decreases with the increase of the path loss 
exponent. This result implies that the energy and spatial 
spectrum efficiency increase with the increase of the path 
loss exponent in Fig. 13 Fig. |7] and Fig. HO] When the path 
loss exponent is fixed, the energy efficiency decreases with 
the increase of the BS density in PVT and grid cellular 
networks. 


V. Conclusions 

To evaluate the spatial spectrum and energy efficiency in 
network level, the Markov chain is first integrated into the 
PVT random cellular networks in this paper. Based on the 
Markov chain based channel access model, spatial spectrum 
and energy efficiency are analyzed for PVT random cellular 
networks. To derive these models, a Markov chain is first 
presented for modeling of wireless channel access in a 
typical PVT cell. Moreover, taking into account the path 
loss and Rayleigh fading effects over wireless channels, the 
outage probability and the blocking probability are derived 
for a typical PVT cell. Furthermore, the spatial spectrum 
and energy efficiency are obtained for PVT random cellular 
networks. Numerical results have shown that the call arrival 
rate in a PVT cell and the BS density of PVT random 
cellular networks have adverse effects on the spatial spec¬ 
trum efficiency of PVT random cellular networks. Moreover, 
the path loss exponent and the SINR threshold have great 
impact on the energy efficiency of PVT cellular networks. 
In the end, our results provide insights into the evaluation 
of spatial spectrum and energy efficiency of PVT random 
cellular networks considering different call arrival rates in 
cellular scenarios. 
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